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All-phononic Digital Transistor on the Basis of Gap-Soliton Dynamics in Anharmonic 
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A conceptual mechanism of amplification of phonons by phonons on the basis of nonlinear band- 
gap transmission (supratransmission) phenomenon is presented. As an example a system of weakly 
coupled chains of anharmonic oscillators is considered. One (source) chain is driven harmonically by 
boundary with a frequency located in the upper band close to the band edge of the ladder system. 
Amplification happens when a second (gate) chain is driven by a small signal in the counter phase 
and with the same frequency as first chain. If the total driving of both chains overcomes the band- 
gap transmission threshold the large amplitude band-gap soliton emerges and amplification scenario 
is realized. The mechanism is interpreted as nonlinear superposition of evanescent and propagating 
nonlinear modes manifesting in a single or double soliton generation working in band-gap or band¬ 
pass regimes, respectively. The results could be straightforwardly generalized for all-optical or 
all-magnonic contexts and has all the promises for logic gate operations. 

PACS numbers: 05.45.-a, 43.25.-hy, 05.45.Yv 


Since the celebrated Fermi-Pasta-Ulam (FPU) first nu¬ 
merical experiment [l[ in 1954, anharmonic oscillator 
chains became a powerful tool in dealing with both fun¬ 
damental aspects of statistical physics [2|, [sf and nonlin¬ 
ear wave phenomena Q and, at the same time, serve 
as the simplest prototypes for extremely complex con¬ 
densed matter systems p, @ and even biophysical pro¬ 
cesses 0 , 0 - In particular, studies on the FPU chains 
together with its further developments, namely nonin- 
tegrable (Klein-Gordon @ and Frenkel-Kont or ova 0) 
and integrable To d a 11 ll c hains, had an impact on the 
discovery of solitons [l2|, ll3| , helped much in understand¬ 
ing of interplay between integrability and chaos 0 , have 
been widely applied for understanding of anomalous ther¬ 
mal conduction and rectification properties in realistic 
physical systems [l5l-[l8| , have applied to describe trans¬ 
port properties in electric transmission lines and even 
in quantum systems, such as Josephson junction paral¬ 
lel arrays and lattices 23 , HBj , and untill now are widely 
used to resolve thermal equipartition issues [IgI . 

Surprisingly, the ladder extension of anharmonic one¬ 
dimensional systems is rarely studied (but see Refs. 

111), although there exists a wide range of appli¬ 
cations for realistic systems, e.g. optical directional cou- 
)lers liS0, weakly coupled classical [l9| or quantum 


21 


in chains, coupled two or multicomponent systems 
[22|, etc. In the present letter we aim to consider 


two weakly coupled FPU chains in order to realize dig¬ 
ital all-phononic amplification of acoustic signals. The 
considered concept of amplification could be straightfor¬ 
wardly extended in case of similar all-optical MM and 
devices. 

ill developments renew the interest 


all-magnonic [33|, 


Phonon laser 35 


in various applications of monochromatic acoustic waves. 
Several ideas have been proposed for phonon diodes 
37 -4^ and all-phononic transistors 13, 45| working on 
magneto-acoustic, nonlinear wave-mixing or mode-mode 


interaction effects. 

In this letter we inrplement a nonlinear band-gap trans¬ 
mission mechanism |46l448| producing gap-solitons in or¬ 
der to achieve digital amplification of weak acoustic sig¬ 
nals. All three ports of the proposed device work on a 
single operational frequency and the schematics is pre¬ 
sented in Fig. la, where the amplifying part of the de¬ 
vice is indicated by a dashed frame. Green and blue color 
chains outside the frame are used for supplying the sig- 



a) 


I Drain 

I k'' K 





FIG. 1: a) The conceptual scheme of phonon transistor with 
indications of source and gate signal forms and supply places. 
Amplified signal is monitored at the Drain port, b) Space- 
time dynamics of gap-soliton creation and propagation in the 
lower chain, where in the range 1 < n < 150 one has a linear 
chain (green balls in the upper graph) and for 150 < n < 250 
soliton propagation occurs in a nonlinear chain (blue balls). 
The upper chain is coupled with the lower one in the range 
150 < n < 200. c) Energies of the input Gate signal (green 
curve) and the amplified gap-soliton output Drain signal (blue 
curve) monitored in the same lower chain 
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nal at the gate and monitoring the output pulse at the 
Drain ports, respectively. The signal is injected from the 
left (gate) linear oscillator chain (green balls). The gate 
signal has a carrier frequency within the band gap of the 
system of the nonlinear oscillator ladder (red and blue 
balls inside the frame) and without source driving can¬ 
not propagate further. At the source input we apply large 
amplitude harmonic driving with the same frequency as 
a gate signal. The amplitude of the source is just below 
the band gap transmission threshold for antisymmetric 
mode of the ladder and together with the gate signal the 
overall amplitude is enough to exceed the threshold and 
a single large amplitude soliton passes the ladder system 
and appears at the drain. While without the gate signal 
the soliton is not produced. Thus it is clear that the soli¬ 
ton amplitude is mostly defined by source driving and the 
digital amplification scenario takes place. Main results of 
the numerical simulations are presented in Fig. lb and 
c). Particularly, graph (b) describes space-time evolu¬ 
tion of the displacements Wn of the lower chain, while in 
graph (c) the energies of signals at the Gate and Drain 
ports are displayed. For clarity we present a movie file 
in Supplemental Material (SM) in order to show the sig¬ 
nal propagation and distribution among upper and lower 
chains. 

For the analytical consideration we examine the ladder 
part of the system (blue and red balls within the frame 
in Fig. la) modeling the system as two weakly coupled 
FPU chains as follows: 

TYlUji = “h Uji—\ ‘lUn) T ^3('^n-|-l "^n) 

+ - Unf + k{Wn - Un) 

mWn = ki{Wn+l + Wn-1 - 2Wn) ^ k2,{Wn+l - Wn)^ 

~\~k2,{Wji—\ 't^n) T k{lLji 'ICn) (1) 

where Un and Wn stand for displacements of the n-th 
oscillators (with mass m) of the upper and lower chains 
respectively; k, k[, /ci, ks and k'^ are linear and nonlinear 
coefficients of stiffness of the springs. Without restricting 
generality we rescale displacement amplitudes and time 
such that the parameters of the upper chain take the unit 
values m = k[ = k';^ = 1. All numerical simulations will 
be done using this scaling and fixing the parameters of 
the lower chain and interchain coupling as follows: ki = 
1.1, ks = 3.5, k = 0.2. We apply dirichlet boundary 
condition at the left end of both chains oscillating the 
balls n = 0 of upper and lower chains with the source 
and gate amplitudes, respectively. 

In the linearized version of o we can readily define 
in the n-th site of the ladder a two component vector 
(un^Wn) and seek for a solution in a form of harmonic 
waves 

iun,Wn) = iR, ( 2 ) 

which gives us two branches of antisymmetric mode 
(neighboring the n-th sites in different chains of the lad¬ 


der oscillate in counter phase) and symmetric mode (the 
n-th sites oscillate in the same phase). Modes with corre¬ 
sponding dispersion relations Ui(p) and ft 2 {p) with max¬ 
imum values at p = TT are displayed in Fig. 2. Those 
modes are characterized by respective components Ri 
and R 2 given by the following formulas 

Rj = k/[k — (Uj)^ + 2(1 — cospj)] (3) 

where j = 1,2. In the case of identical chains in the 
ladder Rj = ±1, i.e. the oscillation amplitudes of in¬ 
terchain neighbor oscillators are the same, while in our 
case of asymmetric ladder the oscillation amplitudes are 
larger in the lower chain (for j = 1) or in the upper 
chain {j = 2). If one drives the ladder with a monochro¬ 
matic frequency Q the excitation wavenumbers pj of the 
respective modes are calculated via the relations: 

U = Ui(p); U = U2(p), (4) 

Depending on the driving frequency (see Fig. 2) the 
following three cases could be realized: 1) for excitation 
frequencies Cl < U2(7r) < Ui(7r) a nonlinear wave en¬ 
ters the system, then separates into two soliton waves; 
2) for excitation frequencies U2(7r) < Cl < Ui(7r) a 
nonlinear wave generates a single soliton associated to 
the antisymmetric mode j = 1; And finally, 3) for the 
Cl > Cli{7r) > Cl2{7r) a nonlinear wave cannot enter the 
system unless the driving amplitude exceeds the band 
gap transmission threshold. 

We start our analysis from considering driving frequen¬ 
cies, which are within the band of both modes (the lowest 
dashed horizontal line in Fig. 2). Then the wavenumbers 
Pj of both modes are real and could be found solving Eqs. 

Then a weakly nonlinear solitonic solution could be 



FIG. 2: Red and blue curves represent linear dispersion 
relations of the antisymmetric Ui(p) and symmetric Cl 2 {p) 
branches, respectively. Dashed horizontal lines represent 
three different cases(D = 1.90, U = 2.10 and Q = 2.18) which 
lead to the different kind of nonlinear wave transport in the 
system. Inset shows the schematics of the system of the FPU 
ladder. 
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FIG. 3: Numerical simulations on the FPU ladder m with the 
boundary driving according to (nnii: the upper chain (upper 
graph) is driven with a frequency U = 1.9 and the lower chain 
(lower graph) is kept pinned. Inset shows how the ultimate 
left n = 0 balls of the ladder are driven in time. 


presented as a modulation of harmonic expression ([2]) of 
the corresponding mode |4^51|: 




where ipj{^,T) is a function of slow variables ^ = e{n — 
Vjt) and T = et^ ^ where Vj = dClj{p)/dp is a group 


p=Pj 

velocity of the respective mode and e is a small expansion 
parameter. At the same time r) obeys the nonlinear 
Schrodinger (NFS) equation (please see for details SM 
file): 




with the following parameters: 

12(1 - cospj)^{ks + Rj) 


A.= 


nj{l + R^A 


“i gp2 


P=Pj 


(6) 


, (7) 


and finally one arrives to the soliton solution of the re¬ 
spective mode as follows: 


= (-Rj, 


1 ) 


Aj cos{ftt — Pju) 
cosh [{n — Vjt)/Aj] 


(8) 


where Aj is a soliton amplitude, while soliton width Aj 
and modified dispersion relation are given by: 


1 2Q'' 1 

Let us note that in the nonlinear case the latter relation 
has to be applied for a computation of soliton carrier 
wavenumber pj. 


FIG. 4: Same as in Fig. 3, but now the upper chain is driven 
by a frequency Q = 2.1 and the lower chain is again pinned. 


In weakly nonlinear limit (small soliton ampli¬ 
tudes Aj ^ 1) and large relative group velocities 

I'^i — V 2 \ /'^i ,2 ^ 1 one can combine the solutions (j8j) ac¬ 
quiring additional phase shift [H^ which could be safely 
neglected in the mentioned limits. By this one is able 
to construct the solution, which describes the initial ex¬ 
citation of the boundary of the solely upper chain. In 
particular, if one takes Ai = A 2 and finds such an ex¬ 
citation frequency that Vi/Ai = 172 /^ 2 , the combination 
— {u^^w‘^) at the origin n = 0 gives 


- {ul,wl) = (i?i -i? 2 , 


0) 


Ai cos(Qt) 
cosh [vit/Ai] 


( 10 ) 


thus driving both chains in time according to the above 
expression one can excite two soliton solution belong¬ 
ing to different branches. That is displayed in Fig. 3, 
driving in numerical simulations the left end of the up¬ 
per chain uq with a frequency ft = 1.90 and amplitude 
A = {Ri — R 2 ) Ai with Ai = 0.025 and calculating Rj 
from Eq. At the same time the lower chain is kept 
pinned at the left boundary {wq = 0) according again to 
the expression m- As seen, the numerical test is just in 
tact with the expectation, as far as according to (|8|) we 
observe different amplitudes for the solitons in the upper 
chain and just the same Ai = 0.025 in the lower one. 

Next we examine one soliton generation driving again 
only upper chain with a frequency lying in the limits 
^2(7r) < ft < f2i(7r), particularly we apply ft = 2.10 in 
numerical simulations (see middle horizontal line in Fig. 
2). In this case antisymmetric mode {j = 1) solution 
could be again presented in solitonic form (jHj), while the 
symmetric mode (j = 2) has no longer a solitonic profile, 
instead it is described by evanescent wave since the corre¬ 
sponding wavenumber p 2 is imaginary number (solution 
of dispersion relation ft = ft 2 {p) has no real roots): 

{ul,wl)={R2, l)B(t)e“lP"l”cos(n2^) (H) 
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FIG. 5: The amplification scenario for FPU ladder. In the 
upper chain (upper panel) we provide the continuous driving 
of the left end with an amplitude just below the threshold 
(HU and the band-gap frequency U = 2.18. While the left 
end of the lower chain (lower panel) is perturbed by a small 
amplitude signal with the same carrier frequency. Inset shows 
how the ultimate left n = 0 balls of the ladder are driven in 
time. 

where can slowly vary in time. This means that we 
observe only one soliton entering the chain. As we try 
to nullify oscillations in the lower chain B{t) should take 
a form of B{t) = Aisech (r^it/Ai) and then the combi¬ 
nation at the origin n = 0 gives the 

same form of the driving as in the previous case ([TQ|) of 
the two soliton generation. The results are displayed in 
Fig. 4, and as seen driving the upper chain with a fre¬ 
quency Q = 2.10 now one monitors the generation of a 
single envelope soliton. 

Finally we consider the case Q = 2.18 (upper dashed 
line in Fig. 2) lying in the band gap of both modes, for 
which only evanescent wave solutions m is realized for 
the modes if the driving amplitude is small. However, if 
the amplitude exceeds some threshold value, a gap soliton 
can be created and propagate along the ladder. For the 
estimation of this threshold value, we assume that the 
upper chain is driven with the amplitude A while the 
lower one is kept pinned. Then, looking at the typical 
solution of such a scenario (nni one can notice that the 
weight of the antisymmetric mode Ai is defined from 
the relation A = Ai (1 — R 1 /R 2 ) and the threshold value 
is calculated from the expression of nonlinear frequency 
shift ([9]): 

A*'* = (1 - y4 [n-fliin)] /Ai. (12) 

Determining gives us an opportunity to realize the 
amplification scenario. For this we create the continuous 
driving in the upper chain with a band-gap frequency 
ft = 2.18 and amplitude just below the threshold, then 
even small counter-phase pulse in the lower chain can 


help to overcome the threshold and provide the necessary 
amplification effect for the weak pulse. For the numerical 
experiment displayed in the Fig. 5 we use a continuous 
driving with the amplitude A = 0.202, while the pulse 
amplitude in the lower chain can be of the order of 0.015. 
As seen such a small pulse is enough to create a gap 
soliton and realize amplification scenario in the oscillator 
ladder. In order to provide a realistic input-output ports 
we have lengthened lower chain adding linear part at the 
left and nonlinear part at the right (see Fig. la) and as 
could be seen the results are in agreement with developed 
analytical scheme. 

This amplification mechanism could be directly ver¬ 
ified using cantilever arrays 0, particularly, one can 
examine two coupled in parallel cantilever arrays and 
use the scenario presented in the Fig. 1. Note that al¬ 
though a cantilever array model include onsite coupling 
terms in contrast to our model of the FPU lattices, the 
consideration of the amplification mechanism will be the 
same, since we consider upper band-gap localized modes 
(staggered excitations) which are similar in both types 
of the anharmonic chains. Using the soft mono-element 
lattice parameters from 0 and taking interchain cou¬ 
pling constant as 20% from onsite coupling coefficient it 
follows that the upper band-gap, appearing due to dis¬ 
creteness, starts at r^i(7r) = 116.8KHz and applying the 
driving in the band-gap with a frequency Q = 117.8KHz 
one can calculate the threshold amplitude for the band- 
gap transmission according to m and it gives the value 
Ath = 0.4/im that is much less than the typical distance 
between cantilever and substrate d = 10/im and excita¬ 
tion can propagate through the array. The width of the 
soliton appearing after the amplification could be calcu¬ 
lated from and gives the value Ai ^ 10 lattice sites. 
Thus it is very clear that the developed amplification 
scheme is very robust with respect to the model choice 
0 - 

Concluding, it should be emphasized that we are us¬ 
ing a single operational frequency and thus the output 
signal could be readily used for the further processing. 
Suggested mechanism could be applied to study amplifi¬ 
cation in quantum systems, e.g. for trapped cold atoms 
in optical lattice ladders. In this case engineering edge 
defect site one can reach the threshold effect by tuning 
the atomic onsite interaction strength. Then above some 
threshold value quantum solitons will be created via reso¬ 
nance process between the edge defect and the quantum 
bound state modes. But this issue needs a further de¬ 
tailed investigation. 
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SUPPLEMENTAL MATERIAL 

We start from the analysis of two weakly coupled FPU 
chains displayed as Eq. (1) in the main text of the 
manuscript: 

TflUji — ‘IVjji) '^n) 

~\~k'^{Un—l '^n) “1“ Un^ 

TflWn — k\{Wfi-^\ “h Wfi—\ “h ^3(^n+l ^n) 

+ /C3(ien-1 - Wn)^ + k{Un - Wn) (13) 


According to a well established procedure |49l-l5lj of 
multi-scaling approach we are seeking a weakly nonlin¬ 
ear solution of ([T3]) in a form of following perturbative 
expansion: 


u = E^° 

a=l 


+ 00 


(14) 


where we define column vector while 

^ and r are slow variables introduced through: ^ = e{n — 
vt) and r = et^; is a soliton group velocity defined 
below and e is a small expansion parameter. 

We go on with equating powers of e substituting ex¬ 
pansion m in set of equations ([T]) . In the linear approx¬ 
imation we have the column vector = 

(/:?(^,r)R and = 0 for \m\ 7 ^ 1; not restricting gener¬ 
ality we can take a space-time independent column vector 
as R = (R, 1), where R is a complex number and 
is a scalar function of slow variables to be determined 
in the next approximations. Then by considering a = 1 
(linear approximation) and the harmonic m = 1 we arrive 
to the equation: 


W*R = 0 


(15) 


where 


W = 


+ 2 (cosp — 1) — k 


k 


k 


0? H- 2ki (cosp — 1) — k 


the solvability of which demands Det(W) = 0, which 
gives us two branches of dispersion relations: 


^ 1,2 = {1 — cosp){lki)k ± 

± y/(l — cosp)^(l — ki^ -|- k‘^ 


(16) 


and two corresponding column vectors Rj = {Rj, 1) with 
Rj expressed with the linear parameters of the problem 
Rj = k/[k — + 2(1 — cosp)], where j = 1, 2. Next we 

introduce a row vector L = (L, 1) through the equation 
L * W = 0, that gives us two row vectors hj. In our case 
the respective components of row hj and column Rj are 
identical Lj = Rj . Thus in linear limit we have following 
matrix relations: 


W (Uj) * Rj = 0, Lj * W (Uj) = 0. 


(17) 


In the following for presentation clarity we omit the in¬ 
dexes j and restore them at the end of the calculations. 

We go on with a second approximation {a = 2) substi¬ 
tuting again m into ([T3j) and considering first harmonic 
m = 1, which leads us to the following equation: 


+ 2i(B - llwI)^R = 0, 


dip 


where 


B = 


sinp 0 


0 ki sinp 
Then multiplying ([T8|) by L one has 

L(B-lia)^R=0. 


(18) 


(19) 


( 20 ) 


In order to identify constant v in the equation above, let 
us take the derivative of ([T^ over p and multiply then 
on the row vector L. One gets: 

L^R = 2Lf^lli-B)R = 0. (21) 

dp \ dp J 

Comparing now o and m we immediately get the 
equality v = dCl/dp, thus the definition for group veloc¬ 
ity, while from (HHI one can solve as follows: 


= -2iW-i(B - llvi)^R. 


( 22 ) 


In the third approximation, equating powers of e for 
a = 3 and first harmonic m = I we have: 

p\-r T ( 2 ) p\ 

WU(3) + 2i(B-llwi)--—+ 2zn^R- (23) 

d^ dr 

- / 92 tj( 1 ) 

-C + 12(1 - cosp)^N|(/5|2(^ ^ g 


where 


C-r'r"'’ ) N=ff I ( 24 ) 

0 V — ki cosp J V ^3 ^ ^ 


Now noting that 


dW ^ la^w ^d‘^n, 


2 (b - S2.i) s C^~- n^I (25) 


2 dp‘^ dp‘^ 


We can further simplify (|M|) multiplying it on L and 
taking into account (j22j) and ([25]): 

T (o—i pgy 

y dp^ ^ dp dp 2 dp^ J d^‘^ 

-l-2iU^LR + 12(1 — cosp)^LN|(/:?|^(/9 = 0 (26) 

dr 
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We can get a final form for ([26]) taking first and second 
derivatives of Eq. (USD over p: 


dW 

dp 


R + W^=0; 

op 


dWdR 

dp‘^ ^ dp dp ^ dp‘^ 


(27) 


In order to estimate the effect we take the following 
approximate values of the problem parameters: m = 
ki = k'l = k2o = = 0.1kg/s‘^^ ks = k'^ = 

10^^kg/s‘^m‘^^ /c4o = k'^Q = 10^kg/s‘^m‘^^ and take weak 
interchain linear coupling coefficient as /c = 0.02kg/s^. 
Then the frequencies of two branches Qi and Q 2 are the 
solutions of matrix dispersion relation 


Solving now d'K/dp from the first equation and substi¬ 
tuting it in the second one and then multiplying it on L 
one gets the following relation: 


dW ^ 1 dW 1 d‘^W 


Det 


+ 2k[ (cosp — 1) — k — k 2 o 
k + 2ki {cosp — 1) 


k 

k - k2o 


= 0 


and now substituting this into the (|26|) and restoring in¬ 
dexes j-s one finally arrives to the Nonlinear Schrodinger 
(NLS) Equation for two nonlinear modes j = 1,2: 

2i-^ + ~ = 0 (29) 

where 


and from the ordinary procedure developed above we 
again get NLS equations with following nonlinear coef¬ 
ficients for the antisymmetric {j = 1) and symmetric 
(j = 2) modes: 


12(1-cospj)‘^{k3 + Rj) 

%(1 + RP ’ 


dp‘^ 


(30) 

p=pj 


and wavenumbers pj are the solutions of respective dis¬ 
persion relations ^2 = AjA^/A. 

We use the same approach considering cantilever ar¬ 
rays. Beginning with modified equations of motion: 


12(1 — cospj)^ (ks + + A:4o + k'/^^R^ 


mUn + k2QUn + - k[{Un+l + Un-1 “ 2 Un) “ 

^3(rf<77,-i-i Uji) k^(^Ufi_i '^n) T k(^Ufi = 0 
TnWji T k20Uri T k/^Qll^ ^l(^n-|-l “t“ ^n—1 ‘2iWfi^ 
fe('?^n+l - - k2,{Wn-l “ Wnf + k{Wn - 'Un) = 0 


Erom this point one can obtain the value of the band-gap 
frequency and the threshold amplitude as shown above. 
In our case ^2i(7r) = 116.8KHz and Ath = 0.4/im with 
the driving frequency Vt = 117.8KHz. 
















